The multi-angle instability in dense neutrino systems 
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We calculate rates of flavor exchange within clouds of neutrinos interacting with each other 
through the standard model coupling, assuming a conventional mass matrix. For cases in which 
there is an angular dependence in the relation among intensity, flavor and spectrum, we find in- 
stabilities in the evolution equations and greatly speeded-up flavor exchange. The instabilities are 
categorized by examining linear perturbations to simple solutions, and their effects are exhibited in 
complete numerical solutions to the system. The application is to the region just under the neutrino 
surfaces in the supernova core. 

PACS numbers: 



1. INTRODUCTION 

There are least two astro-physical contexts in which the 
standard model interaction of neutrinos with neutrinos 
(mediated by Z exchange) could play a important role: 

1) . In the supernova process, either in the region just 
under the neutrinosurface (i.e. the point of last scatter- 
ing, on the average), or in the regions just outside of this 
region. The energy distributions of the various species of 
neutrino in this region are important to the dynamics of 
the explosion, to R process nucleosynthesis downstream, 
and to the observed neutrino pulse on earth, should it 
ever be observed again. These energy distributions may 
be affected greatly by coherent neutrino-neutrino inter- 
actions, even though the effects of ordinary incoherent 
v — v scattering are small. There has been much recent 
work on the combined effects of v-v interactions and os- 
cillations in the region outside the neutrino surfaces for 
neutrinos of different flavors 

2) . In models of the early universe before nucleosynthe- 
sis (but after leptogenesis, etc. and well before v decou- 
pling) in which there are substantial neutrino chemical 
potentials, that is to say, neutrino anti-neutrino imbal- 
ance at much greater level than the electron- positron 
imbalance. If flavors are out of equilibrium as well, then 
v — v interactions become important. When sterile neu- 
trinos are added, these scenarios become arbitrarily com- 
plex. They can easily change the parameters that de- 
termine the He abundance, so that a part of the effort 
expended on their behalf has been to avoid upsetting the 
good fits of the standard theory. The models can be ar- 
ranged to produce models of sterile neutrino dark matter 
in the mass region at least several KeV, coupled with a 
mixing parameter small en oug h not to conflict with cos- 
mic x-ray backgrounds [jj|-|lOl]. 

The present work operates in the context of the super- 
nova problem, although there could well be applications 
to the area b), above, as well. Here we consider some 
effects below the neutrino surface, or, more accurately, 
in the transitional region at about the depth of the av- 
erage last scattering for V{u. T } and vs^x, and below the 
depth of the average last scattering for v e and v e . This 



follows up on previous work , [l2T ] in the same region, 
demonstrating the possibility of flavor instabilities caus- 
ing rapid mixing of the angular distributions and energy 
spectra. In these instabilities, which depend only on for- 
ward v processes, flavor is traded between momentum 
states with no change in the momentum states. The in- 
stabilities arise when neutrino momentum distributions 
are anisotropic. The examples worked out in refs. 
dealing with a region outside of the neutrino-surfaces 
have at their core the type of instability as that described 
11 1, [lij ]. which we characterize as the "multi- angle in- 



stability." 

The present work has three purposes: 

a) To provide a better basis for understanding these 
instabilities, described in the above references simply in 
terms of the plots produced in simulations. In the present 
paper we carry out standard linearized stability analysis 
of the underlying nonlinear systems. 

b) To do a more realistic analysis than was done in ref. 
flU of flavor-scrambling effects in the region just under 
the neutrino surfaces. For example, in the current work 
we show results for a case with 14 bins of v momentum 
space, as compared to 2 bins in ref. 
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c) To point out a second category of possible insta- 
bility, potentially more potent than the one underlying 
the work mentioned above. A feature of the results is 
that small fluctuations of initial distributions away from 
symmetrical cases can have a catalyzing effect, inducing 
complete flavor transfers between up-moving and down- 
moving states. 

Our effects evolve over small distances in homogeneous 
matter, in contrast to MSW resonance flavor trading, 
which depends on passage through more extended re- 
gions with variable electron density. They also act to mix 
up neutrino and anti-neutrino flavors simultaneously, in 
contrast to resonant phenomena, which act on one or the 
other. 

The medium in the region of the neutrinosurfaces has 
density of a few times 10 11 g c -3 at a temperature T = 5— 
7MeV. We shall consider neutrino mixing of two species, 
v e and (following the notation of other authors) v x , where 



is some mixture of v.. 



There are a number of time 
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scales that we can define in this region: 

1. A "very fast" time scale (Gfrij,)" 1 , where is the 
neutrino number density n e ~ T 3 . In our domain this 
time is of order 10~ 3 cm. in units in which c = 1. 

2. The vacuum oscillation time, of order 10 cm. in 
our region, for Sm 2 = .8 x 10~ 4 (eV) 2 . All of the rates of 
change for interesting quantities that we calculate in the 
present paper will be much faster than this rate. 

3. A "fast" rate that is the geometric mean of the 
above two rates. The fact that properties of v systems 
can change at this rate was noted by Kostelecky and 
Samuel [1 31 ] 14\ in their work on self-maintained coher- 
ent oscillations in dense isotropic neutrino gases. More 
recent exploration of this topic can be found in refs. [HI 
and This basic rate recurs in the multi-angle flavor 
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instabilities found in refs. [5( 

4. The rate for scattering from neutrons and protons 
in the medium ~ G^T 3 about ten times as fast as the 
vacuum oscillation rate in the domain of interest. 

The equations for flavor density matrices needed to 
derive our principal results are standard by now 17J . But 
in order to explain some caveats in a later section and 
to set the groundwork for future work we shall briefly 
rederive the standard equations. 

We consider a system of neutrinos interacting with 
each other through the Z mediated neutral current cou- 
plings, and with electrons and positrons in the medium, 
through the W and Z mediated couplings. Interactions of 
j/'s with nucleons in the medium are irrelevant, as far as 
the main developments are concerned, since they are in- 
dependent of v flavor. We also include ordinary neutrino 
oscillations. The effects of ordinary incoherent scattering 
by nucleons and electrons in the medium can be neglected 
over the time scales that we consider. 

We define a e (p) and a x (p), as the respective annihi- 
lation operators for a v & and v x of momentum p, with 
corresponding definitions for the operators that annihi- 
late P's. We introduce the density operators, 



We shall refer to this as the "forward" neutrino-neutrino 
interaction, H v>v , 

H uv (p) = ^ E P , q E {iJ} = e , x [l - cos(tf p , q )] 

x (p»,j(p) - pi,i(p)) (pj,i(q) - pj,M)) 
+(pi,i(p) -Pv(p))(pj,j( t i) . ( 2 ) 

where V is the volume. 

The oscillation (or mass-matrix) term H osc is taken to 
be of the form, 

Hosc(p) = 

Ep |2p| _1 »7[p e) x(p) + Px,e(p) + Pe,x(p) + Px,e(p)] 
+ |2pr 1 ?[p e ,e(p) - Px,x(p) + Pe,e(p) ~ Px,x(p)] 
+G F n e [p eje (p) -/3x,x(p) -Pe,e(p) +Px,x(p)] , (3) 

where the parameters £, r\ determine the standard vac- 
uum neutrino mass parameters according to Sm 2 = 
(£ 2 + V 2 )i tan(0) = T]/£, and the third line in ^ comes 
from the standard model neutrino electron interaction. 

In the main part of what follows, we discuss the time 
evolution of a neutrino system under the influence of the 
Hamiltonian H±, where, 



Hx{p) = H vv {p) + H osc (p) 



(4) 



Since H\{p) in (H|) leaves the momenta of the individual 
neutrinos unchanged, the kinetic energy term (in the ab- 
sence of neutrino mass) is irrelevent to evolution under 
the influence of Hi(p). 

The commutation rules of the density operators are 

[Pi,i(p).Pw(p')] = [Si,lPk,j(p) - <W*,z(p)]<W , 

[Pi,j{p),Pk,l(p')} = [- s i,lPk,j(p) + S o,kPi,l(p)] S v,p' ■ ( 5 ) 

With this apparatus the usual equations for the density 
functions are of the Heisenberg form, 



Pij(p) = «i(p) f a>j(p) 



Pi,j(p) = %(p) t « 1 (p) 



(1) 



where i and j take the values e or x. Note the transposi- 
tion of indices in the definition of the antiparticle density 
operator; the leads to neater formulae below. 

The reason that v-v interactions can change flavor- 
spectrum correlations over a very short time scale is that 
in a pure forward encounter between two neutrinos the 
neutrinos can swap flavors, thus engendering coherent 
(because of forward) effects that are more than phases in 
their wave-functions (because of flavor dependence). We 
start therefore by isolating from the Z exchange coupling 
all of the terms with no momentum exchange, that is, the 
terms that connect a pair of i/'s with momenta p, q to a 
pair with momentum p, q; only the flavors being changed. 



^Pij(p) = [py(p)>#] 



(0) 



We obtain, 



f tPi!j ( P ) = =^ £ q £ fe [pi, k (p)\pkA<ti - pk,M)] 
-pj,fc(p)[p;,fe(q) - Pi,fe(q)] [1 - cos(0 p ,q)] 

+ |p|- 1 [A, P (p)] 4J , (7) 

and, 

|^( P ) = =^ E q E fe [piAp)[Pk,M) - Pk,M] 
-Pi,fe(p)[pi,fe(q) - Pi,fe(q)] [1 - cos(6»p^ q )] 

+ |p|- 1 [A,p(p)] 4J . (8) 
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For brevity we expressed the part coming from the 
commutator with -ff osc in ([3]) terms of a conventional gen- 
erator A(A, £), a 2x2 matrix in the flavor space. 



2. MEAN FIELD EQUATIONS 

In this paper we concentrate on a "mean- field" limit, 
in which the operators (in the occupation number space) , 
Pi j, are replaced by their expectation values in the 
medium, < pij >, after the commutators (in the oc- 
cupation number space) have been performed to give the 
evolution equations (jTj) , dHI) - 

This mean field assumption can be stated as 



(Pi,j(p)Pk,l(p')) = (Pi,j(p))(Pk,l(p')) ■ 



(9) 



In this case we can suppress the indices < pi_j > and 
write the evolution equations ([7]) and ([8]) as matrix equa- 
tions in the flavor indices, generated in canonical fashion 
by effective mode by mode Hamiltonians for the matrices 

Pip), 



Heff(p) 



_ V2Gf 



E q [p(q)-p(q)][i-cos(0 p , q )] 



A . 



(10) 



We emphasize that the formulation based on (flT))) de- 
pends totally on the mean field assumption. Although it 
probably gives the correct answers in the problems that 
we address, we shall return to discuss possible corrections 
to the mean-field approximation. This discussion must 
be based on ([2]) rather than on (|10j). 

For computational purposes we introduce collective 
variables for neutrinos by subdividing the momentum 
space into some small number A^b of regions {a} which 
we designate as beams. For convenience we will do the 
subdivision so that there are the same number of neutri- 
nos + antineutrinos, _/V„, in each beam a, 

N v = E E k»M(P) + MP)] = n v V/N B . (11) 

i P C(q) 

Then we define density matrices for the individual beams, 



pth^ 1 E pmCp), 

P C{q} 



and 



(a) - E Mp), 



PU= N - 



pcM 

so that in each beam, a, we have, 



(12) 



(13) 



(14) 



For calculations we need only two operators in each 
beam for particles and two for antiparticles. We define, 



z = o<°0 



e,e rx,x 



r — o (q) r — 

rx.x 



(15) 



In the two flavor case it has become conventional to write 
the evolution equations for matrices that are hermitean 
3-vectors in the internal flavor space. But it is compu- 
tationally more efficient to use, for each value of a, non- 
hermitian r a 's, introduced above, rather than the two 
hermitian components. 

The effective Hamiltonian is now, 

H = Ea,/3 [ 2 (^ ~ fcV (TP ~ fp) 

+ (z a - z a ){z(3 - zp) x [1 - cos 6 a> p] 
+N V E a [|Par 1 A(r Q +rt+f a + f&) + \p a \~^(z a + z a ) 



(16) 



We have here dropped in ^ the terms that involve only 
J2iPi,i(p) an d J2iPi,i(p)' si nce these quantities are con- 
served under the present dynamics 1 

From ([5] )and (fT5")) the commutation rules of our di- 
mensionless variables are 

[r a , rfa] = N~ x z a S a< p , [r a , zp] = 2N~ 1 r a 5 a ,p 

[r a , fj 3 } = -N~ 1 z a S a< p , [f a , zp] = -2N~ 1 r a S a< p . 

(17) 

The equations for time evolution are given for the flavor 
changing operators by, 

^ = Ep [z<*(rp - fp) - r a (zp - zp) 

x[l - cos6» Q:/3 ] + [r]z a - £r a )/\p a \ + G F n e r a , 



1 dt Ta ~ N rag l^P 

x[l - cos6» Q:/3 ] - (r]z a - £f )/|p Q | + G F n e f a , 
and for the flavor conserving operators by, 

V2Gj-n 
N B 

[1 - cos6> tti/3 ] + 2r]( 



(18) 



i^-z 

L dt^ a 



n b E/j <( f P ~ r p) - r a (r} ~ r*) 



■ r * a )/\v a \ ; 

-^%f^Y:p[rl{rp-rp)-f a {r] 
[1 - cos 6^] - 2r)(f a - r* a )/\p a \ 



(19) 



These terms would enter the dynamics if, for example, we were 
including mixing with sterile neutrinos. 
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Note that in the equations (jT8|) , (|19|) , once we take the 
mean-held limit, the Hermitian adjoints in the original 
operator expressions turn into simple complex conjugates 
of the functions. 

We shall use (|Tg)) . (Til?]) to determine the evolution of 
the system starting from a flavor-diagonal initial state 
at t = 0, that is, from r a (Q) = 0, f a (0) = 0. Consider 
hrst the system with the parameter 77 = 0, and the initial 
values r a = 0, f a = 0. Then the solution has values 
of z a , z a , that are constant in time, with vanishing r a , 
f a . We examine the stability by linearizing around this 
solution. We note that the linear perturbations of the z's 
do not enter (fT8|) at all, so that we have linear equations 
for perturbations Sr a , Sr a . Grouping these perturbations 
into a single column vector 5il a ^ where i = 1 picks out 
Sr a , and i = 2 picks out Sr a , we write the equations in 
matrix form as, 



i—sn = MSQ. 

dt 



(20) 



where M is a 2Ab dimensional square matrix with ele- 
ments M a ,i-pj that we read out from (fT8| as, 

M a ,i- t fi,i = V^GprivNg 1 z a [l - cos6> Qi/3 ] 

7 

M a ,\-j3,2 = - VZGFTlvNg 1 Z a [l - COS 6^] , 

M a , 2;/ 3,i = V^GprivNg 1 z a [l - cos6» Qi/3 ] , 
M a ,i ;) 9,i = -V^Gpn^Ng 1 -z a [l- cos6> Qi/3 ] 

7 

(21) 

The flavor-diagonal initial conditions determine the 2 TVb 
parameters z a , z a in (|21[) . The instabilities that are the 
focus of the present paper all arise from a complex eigen- 
value, Xi of M. There are two categories of instabilities 
that can arise: 

a) Cases in which there is a complex eigenvalue even 
when £ = 0. The growth rate of the mode can be some 
appreciable fraction of Gfii e , the "very fast" rate men- 
tioned in the introduction. This situation only arises 
when the angular distributions in the initial state are 
somewhat complex . 

b) Cases in which the eigenvalues are real when £ = 
but for the case of non-vanishing £ are complex with 
imaginary part of order (Gf^iz/IpI) 1 / 2 , for small £. This 
is the "fast rate" mentioned in the introduction, still 
much faster than the ordinary collision frequency in our 
region of interest. 



3. SIX BEAM EXAMPLE 

We have done a number of simulations with different 
distributions in angle and flavor. Simulations including 
bins of different energies have been done as well, but it 
turns out that including an energy spectrum in the dy- 
namics makes little difference to the results, although en- 
ergy enters in the factor | p> | 1 in the oscillation rate. This 
is basically the synchronization phenomenon discussed in 
refs [HI . Of course, energy enters in a bookkeeping way, 
since the main point will be the equalization of the v x 
spectrum and the v e spectrum through flavor trading. 

Since a flavor independent distribution pij — Si.j is 
completely irrelevent in these equations, we can subtract 
a flavor independent hfi part, replacing the multiplying 
n v in (|18p and (fT9]) by n e g = [n v — ripi), and giving a 
residual distribution that is all (or nearly all) i^'s going 
up and v e 's going down. We begin with one of the simpler 
cases for which we have carried out numerical solutions 
of the evolution equations. 

The physics is that in a critical region the v x , v x have 
an outward bias in their momentum distributions while 
the f e , v e are more istropically distributed. Also there 
is some excess of v e overall , coming from the tail of the 
deleptonization pulse. This excess dies off rather fast 
compared to the time scale of total neutrino energy loss 
from the star. The v x , v x also have a harder energy 
spectrum; as noted above, this latter fact is fairly unim- 
portant as far as the mechanics of flavor exchanges is 
concerned, but of course it is central to the conclusions 
being interesting. 

For the simulation, take six beams, iVg rp = 6, oriented 
in the ±x , ±y , ±z directions in ordinary space, which we 
label as the ±1 ± 2 , ±3 directions. 

We take the following values for the functions in (|18[) . 
(O at t = : 

za(0) - -1 , z 3 (0) = -1 , 
z±i(0) = 9, z±i(0)=0, 
z±2(0) = g, z ±2 (0) = , 
z_ 3 (0) = l + 2 3 , 2_ 3 (0) = 1 . (22) 



Here, the 3 direction is up (i.e., outward) and the choice 
^3«off = Pe[e — Px}x = ~n c g for the density function for 
the upward flux indicates an excess of v x s over z/ e 's. Like- 
wise the part with unity in the down direction reflects the 
excess of v e , v e in the downward direction. A distribu- 
tion that is independent of particle flavor is sterile and 
can be subtracted. We categorize the remainder as: 

1) In the up direction, all v x , D x s. 

2) In the transverse directions, equal numbers of 
Vx 1 &x j v e , but with an excess of v e , measured by the 
parameter g. 

3) In the downward direction, all v ei v e , with a greater 
excess of v e . 



5 



We introduce a dimensionless parameter R to measure 
the ratio of the vacuum oscillation rate to the "very fast" 
neutral current rate, 

R = (<Sm 2 )(GFn (0 )|p|)- x ~ I0~ 9 (23) 

where the estimate is for the conditions that we describe 
above, using the solar neutrino value Sm 2 « 10~ 4 . In our 
numerical studies, we have looked at the region 10 -8 < 
R < 1CT 2 , with fixed, physical sin 2 (26> vac ) = .86. 

In our first example we choose £ to be positive, which 
corresponds to the "inverted hierarchy" . We look at the 
eigenvalues of the matrix M of (|2"T|) in the case g = .2, 
which corresponds to a v e surplus of roughly (y e — D e ) jv e 
of 10% , obtaining, among the twelve eigenvalues, \ one 
complex conjugate pair with imaginary part, 

Im[A] = M[G F R} 1/2 , (24) 

where the exponent 1/2 is essentially exact over a region 
in which R changes by over six orders of magnitude. 

Next we calculate the evolution of the system by solv- 
ing the twelve equations (fT8"|) , (fTTj|) with the initial con- 
ditions given by (|22p . for a range of parameters g (repre- 
senting v e excess), £, with A determined by keeping the 
value sin 2 (29) = fixed. 

In the example shown, we took an electron density cor- 
responding to matter with a density of 5 x 10 11 gc -3 , and 
electron fraction Y e = .4. We take a neutrino density 
n v for each species given by a thermal density at a tem- 
perature of 7MeV; and for definiteness n e s = .2n v , and 
g = .2 as above. 

Again we have calculated the development as a func- 
tion of R over the range TO" 1 > R > 10" 7 . Results are 
shown in fig. 1 for the case R — 1CP 5 . For our range of R 
we find that the first peak corresponding to total flavor 
takeover comes at a mixing time i m i x , 

t-l = G F n eS [R]- 56 , (25) 

where the fit is (within our precision) exact over six 
decades as we change R. We have no insight into this law, 
except insofar as we anticipated something close to R 5 
dependence, based on the eigenvalue dependence noted 
above. Of course we would not have anticipated an exact 
R- 5 law, since the system is completely into a nonlinear 
domain by the time of maximal mixing. 

Now extrapolating to the value R = 10~ 9 appropriate 
to neutrinos with average number density correspond- 
ing to temperature T = 7 MeV, and energies of order 
20 MeV we find that the first turnover of flavor occurs in 
a distance ~ 30 cm. 

Keeping the beam configuration and the initial assign- 
ments fixed as above we have varied the parameters in 
the solutions to see if there are qualitative sensitivities: 

a. Changing the all-over neutrino density, n„ ,(i.e. the 
temperature) or the ratio, n e &/n v , (i.e. the scale of the 
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FIG. 1: Plot of the changing v e , v x asymmetry in the upward 
hemisphere, pi — P2, for the case R = 1CP 4 and inverted 
hierarchy, where {p{i,2}) are the respective densities of the 
two flavors as measured in the units n e s. The initial value, 
pi — p2 = —1 indicates total occupancy of our subset of states 
with v x 's. The unit of time is the inverse of the fast rate 

initial asymmetry) is unimportant. It is true that with 
a smaller asymmetry the flavor exchange will be slower, 
but for the case of our parameter region the exchange 
will still be fast measured in terms of other time scales in 
the problem. In particular the turnover region is small 
compared to the (rather idealized) distance between a 
i/ e -surface and a ^-surface (of last scatterings). 

b. Changing the electron density is also not at all 
interesting; a change by a factor of 2 typically makes a 
change of no more than 10% in the position of the tipping 
point. 

c. For the case of the "normal" hierarchy we get in- 
stability only for very small v e excess. This observation 
agrees completely with the findings of ref . Q . We would 
map out the region of instability in more detail, except 
that the results of the next section, where the initial an- 
gular distributions are taken to be more complex, provide 
a much likelier mechanism for rapid mixing in case of the 
normal hierarchy. 

d. Ordinary scattering is very unlikely in the turnover 
times that correspond to fig. 1. But we could think of 
including it to see where the effect of damping by scat- 
tering leaves the distributions after many, many oscil- 
lations of the type shown in fig. 1. Note that in the 
data shown the time-averaged occupancies retain most of 
the initial asymmetry, in spite of the spikes to reversed 
occupancy. Taking damping from scattering, following 
the prescriptions of McKellar and Thomson [18[ (where, 
in the monoenergetic idealization the effects are just a 
damping term in equations for the off-diagonal parts of 
the density matrix) , we indeed find near zero asymmetry 
and little residual oscillation after only one or two scat- 
tering times. We could not do these calculations using 
the physical values of the parameters, however, because 
this would have required hundreds of the oscillations of 
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which two are shown in fig. 1. However in the next sec- 
tion we find that with more chaotic initial conditions we 
get complete mixing in the time average at shorter times, 
and without scattering. 

4. MORE COMPLEXITY 

We move to better coverage, with 14 rays replacing 
the 6 in the simulations of the preceding section. To the 
earlier configuration of rays we add the eight beams that 
make an angle of 7r/4 with the z axis and are positioned 
in the x, z or y, z planes. For the four of these beams in 
the upward hemisphere, we assign common initial values 
z u ,z u that interpolate between the upward and horizontal 
beams of ((2"2")) . 

z u (0) = -.5 + 5/2 , z u (0) = -1 . (26) 

Similarly for the four new beams in the downward hemi- 
sphere, with common initial values Zj, Zd, where we in- 
terpolate between the completely inward ray in (|22p and 
the horizontal rays to get the common values, 

z d (0) = .5 + 35/2 , z d (0) = -.5 . (27) 

We easily confirm the qualitative conclusions of the 
last system, in this somewhat finer grained simulation of 
the same physics. That is to say, we again find an in- 
stability with growth rate proportional to the geometric 
mean between the "very fast" rate Gp n v and the ordi- 
nary vacuum oscillation parameter, with nearly the same 
coefficient and limitations of parameter domains, and a 
detailed plot recapitulating fig. 1. 

Next we put independent random variations in the 28 
initial values. Unfortunately, a whole landscape of new 
possibilities arise. To quantify this in an example we 
return to the linearized instability condition based on 
the eigenvalues of M of ([21) in our bigger space, but 
now adding independent variations Azi(O), Afj(O), dis- 
tributed randomly on the interval {—.05, .05}. We begin 
by turning off the neutrino mass, Sm 2 — 0. For the case 
of the electron excess parameter, S = .1, we find that 
somewhat over 50 % of the time there is at least one 
complex eigenvalue of M. Note that the hierarchy is not 
an issue at this point since 5m 2 = 0. The v e excess still 
can matter; when we increase the excess v e parameter to 
5 = .2 the probability of a complex eigenvalue decreases 
to around 25%. 

When we proceed to calculate the evolution from the 
full equations (fl"8]) , (fT9|) , we find that in almost every 
case with a complex eigenvalue there is rapid mixing of 
flavors among the rays that go upwards and the rays that 
go downwards, but it is more irregular that the behavior 
shown in fig. 1. We show a typical example in fig. 2 for 
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FIG. 2: For a 14 beam case with the initial conditions as spec- 
ified in (J22J), ((23,(123, with small (10%) initial irregularities 
introduced so that there is a growing mode in the linearized 
equations even in the absence of neutrino mass terms. Plot 
of the changing v e , v x asymmetry in the upward hemisphere, 
Pi — P2 where {/0{i,2}} are the respective densities of the two 
flavors as measured in the units n e g. The initial value — 1 
indicates total occupancy of our subset of states with ^'s. 
The unit of time is the inverse of the fast rate Gftibs- 

the case R — 10~ 4 , where results for both signs of Sm 2 
are plotted. 

We note that the differences between the normal and 
inverted sign are minor. In fig. 3, we show a calculation 
with the same initial conditions but for the value R = 
IQ- 5 . 
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FIG. 3: The same conditions as shown in fig. 2 but with a 
neutrino mass parameter R = 10 -5 , instead of 10 -4 , show- 
ing almost the same behavior for the normal and inverted 
hierarchy cases. 

For even smaller values of R all visible differences be- 
tween normal and inverted hierarchies disappear. For the 
range of values 10~ 6 < R < 10~ 2 we find that (125)) can 
be replaced by, 

t-l = G F n cS [R}- 18 , (28) 

although the fit is no longer perfect. In any case the 
flavor mixing is even faster that that estimated in the 
last section, and it is much more thorough. 
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As an example we take initial conditions defined by be- 
ginning with neutrino number densities characteristic of 
equilibrium at T = 7 MeV, then shift 20% of the down- 
moving v xi v x to up-moving states, and do the reverse for 
the v e , v e distributions. In this case we find total mixing 
in the time corresponding to .2 cm. of neutrino path. In 
general it is clearly not the detailed calculation of evolu- 
tion, as shown in fig. 2, that is important, but rather the 
simple linearized analysis that detects a growing mode. 

The latter depends on the scale of the anisotropies only 
through the fact that the parameter g of ([22]) measures 
the amount of v e relative to the up-down asymmetric 
part of the main v, v distributions. Thus scaled-down 
anisotropies in the presence of a big v e excess lead to 
a smaller chance of instability, according to our earlier 
analysis. 2 We conclude, however, that if the angu- 
lar distributions are anywhere nearly as irregular as the 
usual quantities displayed in the results of simulations of 
the core physics (e.g. Y e ), then there should be many 
opportunities for instability. 



terms enter the analysis in some nastier way? 

That said, and accepting the applicability of the "for- 
ward" Hamiltonian, over time periods for which scatter- 
ing is small, we can ask about the validity of the mean 
field approximation. The essence of the question is cap- 
tured using a truncated form of the neutral current inter- 
action. We return to the effective Hamiltonian (fT6)) and 
discard the oscillation terms and all antiparticle opera- 
tors, and also insert a parameter rj multiplying the term 
in which the operators do not change the flavor , 



H = G^2 [ 2r l r l3 + V z a zp] [1 - COS 9 a ,p] 



(29) 



a,/3 



We take only an up beam A and a down beam B , and 
change to the notation appropriate to a system of N spin 
1/2 particles in group A and N in group B, 



3 ) 



i CA 



i CB 



5. EXAMINATION OF THE THE MEAN FIELD 
APPROXIMATION. 

Although we don't find it obvious from the literature, 
we believe that getting from the underlying field theory to 
the evolution equations for the density matrices involves 
two separate steps that need to be examined separately. 

1. ) Replacing the Hamiltonian by the "forward" Hamil- 
tonian of |(2]). 

2. )Making the "mean field" assumption ((9|). 

The first question is: why do we believe that it is legit- 
imate to ignore almost all of the terms in the neutrino- 
neutrino interaction that correspond (in Born approx.) 
to scattering at a finite angle? Clearly it is because we 
are looking for index of refraction effects, in a sense of or- 
der Gf, rather than scattering, where cross-sections are 
of order G F . 

This oversimplifies; the domains do overlap. For ex- 
ample, in considering neutrinos passing through an in- 
homogenious sea of electrons, density fluctuations on a 
scale much greater than the than the v wavelength, it is 
correct to scatter from a space varying index of refraction 
(coming from the "forward" interaction Hamiltonian) in 
place of explicitly adding up waves coming from the small 
angle scattering on each electron. But are we sure that 
when neutrinos are the target and are themselves sub- 
ject to flavor manipulation from the beam, the index of 
refraction matrix from forward processes embodies all of 
the physics? Or could small angle (rather than forward) 



i<ZA i<ZB 

The (mean field) evolution equations are, 



i—n = G{z 1 r 2 - W\ z i) 
at 

i-^?"2 = G(z 2 ri - r?r 2 Zi) , 

and 



(30) 



(31) 



i-nZi = G{r x r* 2 - rjr 2 ) 
i^Z2 = G{r 2 r\ - r* 2 r x ) . 



(32) 



We take initial 2i(0) = N, ri(0) = 0(spins up) ; 
£2(0) = —N; r 2 (Q) — (spins down). The eigenval- 
ues that determine stability in the linearized system are 
easily read off from (f3Tj) . A = ztGNy^rj 2 — 1. First we 
consider the case 77 = where the eigenvalue shows in- 
stability. We have two groups of spins with no intragroup 
interactions but for each pair, one from one group and 
one from the other, an interaction of the "x-y" model 
form [19( . Now we take the initial values of the spins in 
group #1 all to point in a direction a very small angle a 
from the the z axis, and the spins from group #2 to point 
exactly down. The question is: what is the time scale for 
spin mixing from the above mean field equations? Since 
the growth of the off-diagonal operator is as, 



e Gm a 



(33) 



2 Of course, the v a excess is very time dependent, decreasing 
rapidly during the duration of the neutrino burst. 



The answer is i mix w (GN)- 1 (log a) -1 . 

Next we ask what happens when a — 0, where the 
mean field equations say there is no evolution at all. In 
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ref.[19| we reported complete calculations based on (|29[) 
for cases of N=512 (and fewer) spins in each group. We 
found excellent fits to the form, 

t cv = evolution time to 50% mixed = (GNyHogN (34) 
a result that is backed up by unrigorous analytical work 

Thus although it is technically correct to say that this 
correction to the mean held approximation vanishes in 
the limit of a large number of particles, it is only by a 
logarithm, and logarithms are never huge. That said, the 
comparison of mean-held with an initial tilt a with the 
full a = solution shows that the initial-tilt produces 
more rapid transformations as long as a > N^ 1 . 

This is effectively the case in our present application. 3 
Therefore in our present application we assume that the 
mean field assumptions are justified. We offer a caveat 
here, however; we have not carried out the real 2N spin 
solution when oscillations are present, where, for all we 
know, the outcome could be different than the mean field 
case. It is not guaranteed that the solution will interpo- 
late in a simple way between the two limiting cases that 
we have analyzed. 

We leave this problem for the future, but we can 
sharpen it a little here. Addressing the model of 
in ref. [ljj, we also reported results for the case T] = 1, 
where the eigenvalues A are real. Here we found was 
no speed-up of the form (f34|) . Instead we found, t cv ~ 
G-ijV-V 2 , slower by a factor of y^/logN. This can 
be characterized as the perturbative time, since for small 
times, the transition probability for a particular spin is 



prob ~ G 2 Nt 2 



(35) 



Thus our conjecture is that when the linear response cal- 
culation gives an exponentially growing mode, then the 
real solution gives speed-up as defined in (j34|) . In con- 
trast, when the eigenvalues, A, are real, we get the per- 
turbative rate. 

More evidence for the latter is found in Friedland and 



Lunardini's analytic solution [2lJ of a model of the latter 
category (real eigenvalues for the perturbations in the 
mean field equations). Their model was basically that of 
(|29|) with r) = 1 together with the assumption of isotropy, 

[l-COS0 a ,/j]-»l. 4 



3 We assume that the effective value of N would be the number 
of neutrinos in a volume of dimension of some neutrino coher- 
ence length, not that we know exactly how the latter should be 
denned. 

4 The same solution technique can be used to solve the case of 
r\ = 1 with two opposed beams that we discussed above, and also 
the complementary case of r\ = and isotropy, in which case the 
mean field eigenvalues are A = {0 , zi(0) + 22(0)}, real for any 
initial configuration. The present remarks apply to the solutions 
of similar models in rcfs 22] as well. 



6. DISCUSSION 

We have made a case for rapid spectrum (and angle) 
mixing for systems that have two kinds of neutrinos and 
antineutrinos. We fully expect the case with all three 
flavors, and more mixing parameters, to mix rapidly in 
everydirection; a primitive example was worked out in 



ref. [11 1 



To compare with other authors' results in related cal- 
culations, first we emphasize that most published exam- 
ples of models with flavor exchange depend on MSW 
transitions over a distance in which the electron density 
is changing significantly. By contrast, our instabilities 
and rapid flavor exchange do not depend on changing 
the electron density and in any event can occur in regions 
with parameters that are far from resonance parameters, 
as in the region just under the neutrino-surface that we 
focused on here. The transition distances that we find 
in our calculations, of order 1 cm., are also much shorter 
than the resonance regions that appear at smaller elec- 
tron densities. 

The basic mechanism at work in our calculation is sim- 
ilar to that found in refs. 0-0 > which consider rapid 
flavor exchanges just downstream from the neutrino- 
surface. Though it is not so far from the region that we 
considered, the neutrino angular distributions are very 
different in the two regions. Outside of the neutrino sur- 
face all of the flow is outwards and the relative angles 
between neutrino trajectories are diminishing rapidly as 
one moves outwards. Since these angles are the key to 
all v — v effects (as they are in our region as well) the 
results cannot be compared in detail with ours. 

But we can compare general features of the results. 
The time scale for the rapid flavor exchanges found in 
refs. [|| and is the "fast" rate as defined in our intro- 
duction, the geometric mean of the neutral current rate 
and the vacuum oscillation rate. As in our case the phe- 
nomenon does not depend on the space variation of the 
electron density. Except for the geometry, the physics 
appears to be the same as in our "six beam" example. 
We confirm the conclusions of ref. @ as to the critical 
roles (and linking) of the v e excess and the normal vs. 
inverted hierarchy choice. In the case of normal hier- 
archy, the exchange occurs only for very small values of 
the electron excess. The case of inverted hierarchy allows 
considerable electron excess however. 

Our results in sec. 4, where we considered a less sym- 
metrical angular distribution, are quite different. When 
the linearized stability analysis in the absence of neutrino 
mass terms shows a growing mode, then the solution for 
the complete problem shows large mixing for both the 
normal and inverted hierarchy cases, in the presence of 
moderate v e excess. The evolution rate is also an order 
of magnitude or more than in the more symmetric "six 
beam" example. 
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Although the authors of ref. [B} did not find such be- 
havior in the region that they treated, it is possible that 
if perfect cylindrical symmetry in the initial distributions 
were given up, then these super-fast turnovers would be 
produced for both the normal and inverted hierarchy. 

We could ask how much it matters whether or not the 
emerging energy spectra are partially or fully homoge- 
nized. There is strong evidence in the literature [23rj-[25l| 
that it would matter significantly for the R process nu- 
cleosynthesis yields, though we believe that there is no 
single graph showing the effects with precision, as there 
is other complex physics involved. An interesting ques- 
tion is that of whether or not it matters to the actual 
explosion dynamics. The neutrino heating in the region 
above the neutrino-surfaces is augmented appreciably in 
a scenario in which v e 's get boosted in average energy by 
trading spectra with v x 's. 

Of course, feeding our detailed considerations into a 
big code for the supernova is out of the question. The 
hydrodynamics plus neutrino transport problem has not 
been formulated in a way that admits collective effects 
of the kind we have discussed here, so far as we know; 
and there seems little possibility that it could be in a way 
that allowed whole-star computations in any reasonable 
reasonable amount of computer time. 

But it would be easy for the numerical simulators to 
check the limiting case of our suggestion, total instanta- 
neous flavor homogenization everywhere within the outer 
neutrino-surface, putting it by hand into the codes. 

The third consequence, and most obvious observable 
consequence of our considerations, is the effect on the 
neutrino pulse signals from a nearby supernova. Much 
has been written on this subject, but the chance of ob- 
servational data within the near future is very small, and 
loose ends can be cleared up after the fact, since the de- 
sign of observing apparatus appears not to depend on the 
details of the predictions. 

This work was supported in part by NSF grant PHY- 
0455918. 



[S 
[9. 
[10 



[11 
[12 

[is; 

[14 

[is; 
[ie; 
[17; 
[is; 

[19 
[20 
[21 
[22 

[23 
[24 
[25 



C. T. Kishimoto and G. M. Fuller, larXiv:0802.3377l 
X. Shi, G. M. Fuller, Phys.Rev.Lett. 82, 2832 (1999) 
K. Abazajian, N. F. Bell, G. M. Fuller, Y. 
Y. Y. Wong,Phys.Re v. D72, 063004 (2005) , 
|arXiv:astroph/0410175"1 

R. F. Sawyer, Phy s.Rev. D72, 045003 (2005), 

|arXiv:hep-ph/0503013| 

R. F Sawyer, |arXiv:hep-ph/0408265l 

A. Kostelecky and S. Samuel, Phys.Rev. D52 , 621 
(1995) 

S. Samuel, Phys.Rev. D53, 5382 (1996) 

G.G. Raffelt, AYu. Smirnov, Phys. Rev.D76, 081301 

(2007), arXiv: 0705.1830 

S.Hannestad, G.G.Raffelt, G.Sigl, Y.YY.Wong, 
Phys.Rev. D74, 105010 (2006), |arXiv:astro-ph/0608695| 
S. Pas tor and G. Raffelt, Ph ys.Rev.Lett. 89, 191101 
(2002), |arXiv:astro-ph/0207281| 

B. H.J. McKellar and M.J. Thomson, Phys.Rev. D49, 
2710 (1994) 

R. F. Sawyer, Phys.Rev. 

(2004) |arXiv:quant-ph/0312217l 
R. F. Sawyer, Phys.Rev.Lett. 93, 



A70,022308 



133601 (2004), 
JHEP 0310, 043 (2003), 



|arXiv:hep -ph/0404247 
A. Friedland and C. Lunardini, 
|arXiv:hep-ph/0307140| 

A. Friedland, B~ H. J. McKellar and I. Oku- 
niewicz, Journal-ref : Phys.Rev. D73, 093002 (2006), 
|arXiv:hep-ph7 0602016 

R. D. Hoffman , S. E. Woosley and Y.-Z. Q ian, Astro- 
phys.J. 482, 951 (1997), arXiv: |astro-ph/9611097| 
Y.-Z. Qian and S. E. Woosley, Astrophys.J. 471, 331 
(1996), |arXiv:astro-ph/9611094| 

B. S. Meyer and J. S. Brown, Astrophys. J, Suppl 12, 
199 (1997) 



[1] G. M. Fuller, Y.-Z. Qian Phys.Rev. D73, 023004 (2006), 

|arXiv:astro-ph/0505240"1 
[2] H. Duan, G. M. Full er and Yong-Zhong Qian Phys.Rev. 

D74, 123004 (2006), [arXTv:astro-ph/0511275| 
[3] H. Duan, G. M. Fulle r, J Carlson, Y.-Z~Qlan Phys.Rev. 

D74, 105014(2006), [arXiv:astro-ph/0606616| 
[4] H. Duan, G. M. Fuller, J. Carlson, Y.- 

Z. Qian, Phys.Rev. Lett. 97, 241101 (2006), 

arXiv:astro- ph/0608050| 

[5] GTGRaffelt, G.Sigl , |arXiv:hep-ph/0701182| 

[6] A. Esteban-Pretel, S. Pastor, R. Tomas, G. Raffelt and 

G. Sigl, Phys.Rev. D76, 125018 (2007), arXiv: astro- 

ph/0706.2498 

[7] H. Duan, G. M. Fuller, J. Carlson and Y.-Z. Qian, Phys. 
Rev. D75, 125005 (2007), [arXTv:astro-ph/0703776| 



